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ABSTRACT 

We present the first numerical simulations in coupled dark energy cosmologies 
with high enough resolution to investigate the effects of the coupling on galactic and 
sub-galactic scales. We choose two constant couplings and a time-varying coupling 
function and we run simulations of three Milky-Way-size halos (^IO^^Mq), a lower 
mass halo (GxIO^^Mq) and a dwarf galaxy halo (SxIO^Mq). We resolve each halo 
with several millions dark matter particles. On all scales the coupling causes lower 
halo concentrations and a reduced number of substructures with respect to ACDM. 
We show that the reduced concentrations are not due to different formation times, 
but they are related to the extra terms that appear in the equations describing the 
gravitational dynamics. On the scale of the Milky Way satellites, we show that the 
lower concentrations can help in reconciling observed and simulated rotation curves, 
but the coupling values necessary to have a significant difference from ACDM are out¬ 
side the current observational constraints. On the other hand, if other modifications to 
the standard model allowing a higher coupling (e.g. massive neutrinos) are considered, 
coupled dark energy can become an interesting scenario to alleviate the small-scale 
issues of the ACDM model. 

Key words: cosmology: dark energy, dark matter - galaxies: haloes, formation, 
evolution - methods: numerical 


1 INTRODUCTION 

Since the discovery of the accelerated expansion of the 
universe (Riess et al. 1998; Perlmutter et al. 1999), a Cos¬ 
mological Constant A has been the most widely accepted ex¬ 
planation for the required negative pressure. Together with 
cold dark matter, today the dark sector is accounting for 
about 95% of the total energy density (Planck Collabora¬ 
tion et al. 2015) and builds the foundations for the so-called 
A Cold Dark Matter (ACDM) model. 

Despite the highly successful inflationary ACDM 
paradigm, the fundamental problems associated with the 
introduction of a cosmological constant, namely fine-tuning 
and coincidence problems (Weinberg 1989), have served as 


* penzo@mpia.de 


motivations for alternative descriptions of the dark sec¬ 
tor. Introducing a time evolving scalar held (dark energy) 
responsible for the negative pressure is the approach of 
quintessence models (Wetterich 1988; Peebles & Ratra 1988) 
and has been one of the most popular generalizations for 
the cosmological constant in the last decade. Furthermore, 
given the currently still unknown nature of the dark sector, 
the possibility of a non-null coupling between dark matter 
and dark energy has been considered (Wetterich 1995; An¬ 
derson & Carroll 1998; Amendola 2000; Billyard & Coley 
2000; Zimdahl et al. 2001; Amendola & Tocchini-Valentini 
2001; Farrar & Peebles 2004; Gromov et al. 2004). Given 
that in these models dark matter and dark energy density 
evolutions are strongly coupled, this would in turn alleviate 
the coincidence problem (Mangano et al. 2003; Matarrese 
et al. 2003). The effects of such interaction might be seen 
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on the Cosmic Microwave Background (CMB), on super¬ 
novae and on the growth of structures, as pointed out by 
Matarrese et ah (2003); Amendola (2004); Amendola et ah 
(2004); Koivisto (2005); Guo et ah (2007) and many oth¬ 
ers. Structure formation has been as well investigated via 
numerical simulations by Maccio et al. (2004); Baldi et al. 
(2010); Li & Barrow (2011); Carlesi et al. (2014) and their 
follow up works, where the statistical distribution of struc¬ 
tures has been studied. Both Baldi et al. (2010) and Car¬ 
lesi et al. (2014) found that, when introducing a coupling 
between dark energy and dark matter, halo concentrations 
decrease. 

In this work, we run the first high resolution simula¬ 
tions on galactic scales in coupled dark energy cosmology. 
Our aim is to obtain high enough resolutions to investigate 
the properties of the dark matter distribution at sub-galactic 
scales, mass scales at which the effects of the coupling have 
not yet been studied. The subhalos that we are interested in 
will in turn be the hosts of dwarf galaxies and their prop¬ 
erties can be compared with observations of satellite dwarf 
galaxies of both Milky Way and Andromeda. In fact, despite 
ACDM predictions on large scales being in very good agree¬ 
ment with galaxy clustering surveys (Jones et al. 2009; Alam 
et al. 2015), on galactic scales challenges between ACDM 
predictions and observations have appeared. 

Firstly, the missing satellites problem, i.e. overabun¬ 
dance of substructures in ACDM Milky-Way size halo sim¬ 
ulations when compared to observations of the Milky Way 
dwarf galaxies (Klypin et al. 1999; Moore et al. 1999). On 
the other hand, as showed in Madau et al. (2008) and Maccio 
et al. (2010), accounting for the baryonic physics drasti¬ 
cally reduces the number of visible satellites. Secondly, the 
eore/eusp problem, namely the inconsistency between the 
constant density cores estimated from observations and the 
cuspy inner density profiles found in ACDM simulations. 
See Flores & Primack (1994); Moore (1994); Diemand et al. 
(2005); Gentile et al. (2009); Walker & Penarrubia (2011); 
Agnello & Evans (2012); Salucci et al. (2012), but also van 
den Bosch & Swaters (2001); Swaters et al. (2003); Simon 
et al. (2005). While this inconsistency can be attributed to 
baryonic feedback processes (e.g. Governato et al. 2012; Di 
Cintio et al. 2014; Ohorbe et al. 2015), for the case of Milky 
Way satellites the baryonic explanation is not straightfor¬ 
ward since these objects can be almost completely dark mat¬ 
ter dominated. Baldi et al. (2010) and Carlesi et al. (2014) 
showed that for halos with M > 10^^M© the coupling be¬ 
tween dark matter and dark energy produces density profiles 
that are less cuspy in the inner density regions, which can 
help alleviating the core/cusp problem. The aim of this work 
is to investigate whether this effect persists at much lower 
masses. Moreover, concentrations of the most massive sub¬ 
halos orbiting around a ACDM Milky-Way size halo seem 
to be too high to be hosting the brightest dwarf galaxies 
observed. This translates into a prediction from ACDM nu¬ 
merical simulations for the existence of massive dark matter 
subhalos that seem to have failed at forming stars, and is 
known as the too big to fail problem (Boylan-Kolchin et al. 
2011; Lovell et al. 2012; Rashkov et al. 2012; Tollerud et al. 
2012). 

Whether these issues bring serious challenges for the 
ACDM model or whether they can entirely be treated by 
invoking baryonic physics is currently under debate. With 


this work we aim at studing the properties of halos and their 
sub-structures to determine whether coupled dark energy 
cosmologies can alleviate the aforementioned issues. In Sec¬ 
tion 2 we summarize the theoretical model behind coupled 
dark energy and we specify our choices of coupling functions. 
In Section 3 we described the numerical methods used to 
produce initial conditions and the N-body codes to run the 
simulations. In Section 4 we show our simulations results, 
for both halos and subhalos. Finally, in Section 5 we present 
our conclusions. 


2 COSMOLOGICAL MODELS 

We present a study focused on understanding the non¬ 
linear effects of coupled dark energy models on galactic 
scales. The models that we consider allow for an interac¬ 
tion between dark matter and dark energy (Amendola 2000; 
Billyard & Coley 2000; Zimdahl et al. 2001; Gromov et al. 
2004; Maccio et al. 2004; Baldi et al. 2010) and obey the 
following Lagrangian: 

^ + CkinW, (1) 

where the mass m(0) of the dark matter field '0 is a function 
of the dark energy scalar field 0, and the Ckin[ip] term in¬ 
cludes the kinetic part of the dark matter Lagrangian. The 
choice of m(0) specifies the coupling and in our work we use: 

— moe ^pi , (2) 

where mo is the mass at z = 0, Mpi = 1/ VStvG is the Planck 
mass, with G being the Newton’s constant, and /3((/)) is the 
coupling function. The respective continuity equations for 
cold dark matter and dark energy are: 

pc + 3iLpc = -I3((l))(j)pc, (3) 

p(j) + 3Hp(f) = +/3(0)0pc,, 

where pc is the cold dark matter density and is the dark 
energy density, which is pcj^ = + 1/(0), and if = a/a is the 

Hubble parameter. Our choice for the self-interacting dark 
energy potential is 1/(0) oc e~^^^ with a = 0.08. 

The evolution of cold dark matter density perturbations 
is regulated by the following equation: 

5c + {2H - /?^)5c - [(1 + 2l3^)ncSc + = 0, (4) 

where Qc and are respectively the density parameters 
Qi = pifpcrit for cold dark matter and baryons, and pcrit = 
SttG/SH^. Two extra terms appear in Eq. 4 compared to the 
ACDM case: a friction term —^(j)6c and the factor (1 + 20^) 
responsible for the enhancement of the gravitational force 
acting on cold dark matter particles, which is known as “fifth 
force”. As pointed out in Baldi (2011), in the linear regime 
both these extra terms produce an acceleration of growth of 
cold dark matter density perturbations. On the other hand, 
when considering the non-linear effects, the friction term is 
responsible for lowering the concentration of dark matter 
halos. 

The appearance of extra terms becomes clear when cal¬ 
culating the acceleration felt by the z-th dark matter particle 
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a = l/(l+«) 


Figure 1. Linear growth factor evolutions for all cosmologies 
normalized to today’s values divided by the ACDM evolution. 


in a coupled dark energy cosmology for the limit of a light 
scalar field (see Baldi et ah 2010 for calculation): 

Vi — l3{(j))^Vi + G[1 + 2^(0)^] ^ I - 7/ • (5) 

The term accelerates dark matter particles in the 

direction of their motion and thus lowers halo concentra¬ 
tions. 

Based on Baldi et al. (2010) and Baldi (2011), we chose 
three coupling scenarios. EXP003 and EXPOOSeS are obser- 
vationally viable models (see Pettorino et al. 2012 for CMB 
constraints on the coupling value), respectively with a con¬ 
stant coupling /3 = 0.15 and a coupling /3(0) varying with 
redshift (for more details on these models we refer to Baldi 
2012). We also explore an extreme constant coupling case 
EXP006 with {3 — 0.3 (about 6cr outside present observa¬ 
tional limits, Pettorino et al. 2012) to better understand its 
implications. 


3 NUMERICAL METHODS 

3.1 Initial Conditions and Coupled Dark Energy 

As in Penzo et al. (2014), we used grafic-de, 
an extension of the initial condition generator GRAFIC-2 
(Bertschinger 2001) such that initial conditions for a generic 
cosmological model can be produced once the evolution of 
the cosmological parameters are given as an input. GRAFIG- 
DE requires transfer functions, evolution of the density pa¬ 
rameters Qi, linear growth factor D+ and fn, the logarith¬ 
mic derivative of the growth factor with respect to the scale 
factor. As the original code, grafig-de is able to generate 
multi mass initial conditions from a cosmological box. In 
Eig. 1 we show the evolution of the linear growth factor D+ 
for all four cosmological models. The transfer functions for 
ACDM have been produced using GAME (Lewis & Bridle 
2002), while the transfer functions for the coupled dark en¬ 
ergy models have been produced by weighting the ACDM 
transfer functions with the D+ of the coupled model. Ah ini¬ 
tial conditions were created using the same random seeds, 
in order to be able to identify structures among the models. 



Figure 2. Ratio between the mass functions for the coupled dark 
energy cosmologies and the one from ACDM for redshifts 2 : = 
2,1, 0 for the cosmological boxes of size 80 Mpc/h. 


All the models share the same cosmological parameters and 
have at z = 0: = 0.0458, ^Idmq = 0.229, Ho = 70.2 km 

s“^ Mpc“^, (78 = 0.816, ris — 0.968, where these parameters 
are density parameters for baryons and dark matter, Hub¬ 
ble constant, root mean square of the fluctuation amplitudes 
and primeval spectra index. 


3.2 N-body Simulations 

We hrst generate two sets of uniform particle distribu¬ 
tions, a 80 Mpc/h box and a 12 Mpc/h box, both with 350^ 
particles. The initial conditions were evolved with the code 
GADGET- 2 (Springel 2005), which includes the coupled dark 
energy implementation introduced in Baldi et al. (2010). In 
Fig. 2 we show the ratio between the mass functions for the 
coupled dark energy cosmologies with the one from ACDM 
for redshifts z = 2,1, 0 for our 80 Mpc/h boxes. 

We chose four dark matter halos in the ACDM 80 
Mpc/h box and one dwarf halo in the ACDM 12 Mpc/h box, 
and looked for their corresponding realizations in the cou¬ 
pled dark energy simulations. Note that we used the same 
random seed for all initial conditions to be able to follow 
the same halos in ah cosmological boxes. Our halos have 
been chosen so that no other halos with comparable masses 
were found within four of their virial radii. We then re-ran 
the cosmological boxes with increased resolution in a La- 
grangian volume that includes ah particles that at 2 ; = 0 
were found in three virial radii of each selected halo. 

Our hnal sample is composed of three Milky-Way-size 
halos (haloo:, haloy^ and haloy), a 6x lO^M© halo (halo(5) 
and a dwarf halo (haloe). For more details on the halos 
properties at z = 0 see Table 1. For the halo identihca- 
tion we used the code Amiga Halo Finder (Knohmann & 
Knebe 2009). The softening lengths are chosen to be 1/40 of 
the intra-particle distance in the low resolution simulation 
divided by the rehnement factor RF; RF = 15 for haloo:, 
halo/3 and haloe, RF = 24 for haloy and halo(5. Precisely, 
the softening lengths are 0.54 kpc for haloo: and halo/3, 0.34 
kpc for haloy and halo(5, 0.081 kpc for haloe. The parti¬ 
cle masses at z = 0 in the high resolution volumes are 
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Table 1. Physical properties of the five halos in all cosmologies, 
ACDM, EXP003, EXP008e3 and EXP006. We show mass at R 200 , 
R 200 , concentrations and number of particles within R200- 



M 200 
[Mq 1 

R 200 

[kpc] 

c = R 20 o/rs 

N 200 

haloo 

ACDM 

2.6x10^2 

284 

11.8 

6.8x10® 

EXP003 

2.5x1012 

281 

9.1 

6.6x10® 

EXP008e3 

2.6x10^2 

282 

10.2 

6.7x10® 

EXP006 

2.1x1012 

265 

4.6 

5.5x10® 

halo/3 

ACDM 

2.5x1012 

278 

10.7 

6.3x10® 

EXP003 

2.2x1012 

267 

8.0 

5.7x10® 

EXP008e3 

2.2x1012 

268 

8.7 

5.8x10® 

EXP006 

1.7x1012 

246 

4.3 

4.5x10® 

haloy 

ACDM 

9.7x1011 

204 

10.8 

1.0x10'^ 

EXP003 

9.3xl0ii 

201 

8.6 

9.9x10® 

EXP008e3 

9.5x1011 

203 

9.6 

1.0x10'^ 

EXP006 

7.6x1011 

188 

3.2 

8.1x10® 

halo5 

ACDM 

6.4xl0ii 

177 

13.3 

6.8x10® 

EXP003 

5.6x1011 

170 

9.7 

6.0x10® 

EXP008e3 

5.9xl0ii 

172 

11.0 

6.3x10® 

EXP006 

5.3x1011 

166 

4.7 

5.6x10® 

haloe 

ACDM 

4x109 

33 

15.3 

3.1x10® 

EXP006 

3x10® 

30 

6.6 

2.4x10® 


S.SxlO^M© for haloa and halo/3, O.dxlO^M© for haloy and 
halo(5, 1.3 xlO^M© for haloe. In Fig. 3 we show the pro¬ 
jected density maps of the four most massive halos for each 
cosmological model. We will discuss the dwarf halo in Sec¬ 
tion 4.3. For the density maps and throughout the paper 
we chose to calculate halo properties using R 200 , radius at 
which the density is 200 times the critical density Pc, with 
p. = 3//V(8.G) 


4 RESULTS 

4.1 Host Halos Properties 

We first analyze the properties of the four most massive ha¬ 
los, haloo, halo/3, haloy and halo(5. We show concentrations, 
density profiles, rotation curves, evolution of the scale radius 
and accretion histories. 


4 . 1.1 Concentrations and Density Profiles 

By introducing a coupling between dark matter and 
dark energy, halo concentrations decrease. This was also 
shown in Baldi et al. 2010, Li & Barrow 2011, and Carlesi 
et al. 2014 for halos with masses M > lO^^M©. In this work 
we investigate mass scales M < lO^^M©. The resolution 
that we are able to reach is higher thanks to the multi-mass 
technique. In Table 1 we show the concentration values for 
each halo, for which we use the definition 

c = R2oo/rs, ( 6 ) 



radius [kpc] radius [kpc] 


Figure 5. Rotation curves for haloo, halo/3, haloy and halo(5 at 
z = 0, each for ACDM , EXP003, EXP008e3 and EXP006. The 
inner radius is equal to three times the softening length, while the 
outer radius is four times R 200 of each halo. The vertical dashed 
lines mark R 200 for each halo in each cosmology. 


where R 200 is the radius at which the density equals 200 
times the critical density and r^ is the scale radius in the 
Navarro Frenk and White (NFW) profile (Navarro et al. 
1997). We computed r^ via a minimization procedure us¬ 
ing the Levenberg & Marquart method as in Maccio et al. 
(2008). In agreement with literature, we find that halos 
which lived in a coupled dark energy cosmology have lower 
concentrations. Fig. 4 shows the density profiles for haloa, 
halo/3, haloy and halo(5. The behavior of the profiles as a 
function of cosmology is maintained for all four halos, with 
a significant flattening of the inner part of the profiles only 
for the extreme coupled cosmology EXP006, while differ¬ 
ences are less evident in the EXP003 and EXP008e3 ha¬ 
los. Interestingly, the EXP006 realization of haloy (M = 
T.GxlO^^M©) produces a much flatter halo profile, with 
slope (a=-0.8, which falls out of NEW parametrization. On 
the other hand, all other profiles of haloa, halo/3, haloy and 
halo(5 in all cosmologies are well described by the NEW pro¬ 
file. 

Additionally, in Eig. 5 we show the rotation curves at 
z = 0 for the four halos. Eor models within the observa¬ 
tional constraints the rotation curves are not significantly 
affected. The only case in which we observe a considerable 
flattening is the extreme model EXP006, for all four cases. 
This is in agreement with Penzo et al. (2014), where we find 
that differences in rotation curves among models within ob¬ 
servational constraints for dynamical dark energy are not 
significant in the dark matter only case. On the contrary, in 
hydrodynamical simulations we find observable differences 
in rotation curves due to the effects of baryons which en¬ 
hance the variations in the dark matter accretion. We expect 
the same enhancement also in coupled dark energy models 
once hydrodynamics is taken into account. This aspect will 
be explored in a future work. 
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Figure 3. Projected density maps of our sample at 2 : = 0. From first row to last we are showing haloa, halo/3, halo7,halo(5; from first 
column to last we are showing ACDM, EXP003, EXPOOSeS and EXP006 cosmologies. All images are spheres of radius R 20 O 5 radius at 
which the density is equal to 200 times the critical density. 


4 . 1.2 NFW Scale Radius Evolution 

In the Section 4.1.1 we have showed that halos that 
form in a coupled dark energy cosmology with a high value 
for the coupling constant have concentrations that are sig¬ 
nificantly lower at z = 0. Given that almost all halos are 
well described by a NFW density profile, it means that their 
NFW scale radii are much larger than the scale radii of 
the corresponding ACDM realizations. In Fig. 6 we show 
the behavior of the scale radius Vs as a function of redshift 
for halo/3 in all four cosmologies; the other Milky-Way size 


halos have similar behaviors. Compared to the ACDM case, 
halos which live in coupled dark energy cosmologies show a 
larger scale radius at all redshifts. 

4 . 1.3 Main Halos Accretion Histories 

In order to investigate the origin of the different concen¬ 
trations, especially in the EXP006 cosmology, in this Section 
we concentrate on the halo formation times. Firstly, in Fig. 7 
we show the accretion histories, namely the evolution of the 
mass enclosed in R 200 normalized to the value of the mass 
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Figure 4. Density profiles for haloo:, halo/3, haloy and halo(5 at 2 ; = 0, each for ACDM (black), EXP003 (cyan), EXP008e3 (blue) and 
EXP006 (red). The inner radius is equal to three times the softening length, while the outer radius is four times R 200 of each halo. The 
vertical dashed lines mark R 200 for each halo in each cosmology. 



Figure 6. Scale radius obtained by fitting an NEW density profile 
using the Levenberg & Marquart method for halo/3 as a function 
of redshift. 

at z = 0 as a function of expansion factor a = 1/(1 + z). 
Halos growing in ACDM, EXP003 and EXP008e3 cosmolo¬ 
gies show similar accretion histories especially for haloa and 
haloy; while halos living in the EXP006 cosmology accrete 
their mass earlier on compared to their ACDM realizations. 
Among the three halos, coupled cosmologies runs show un¬ 
expected drops in the accretion histories. These would be 
unusual in a ACDM scenario since halo total masses do not 
decrease unless it is a temporary effect of a merger (see for 
instance haloa and halo/3 around a = 0.3). On the other 



Figure 7. Evolution of the mass enclosed in R 200 normalized to 
the mass at a = 1 as a function of the scale factor for haloo, 
halo^, haloy in all four cosmologies. 


hand, in the case of coupled cosmologies the friction term 
in Eq. 4 is responsible for injecting kinetic energy into the 
system, which may cause some particles to become gravita¬ 
tionally unbound. 

In order to estimate the time of formation for each halo, 
we followed the approach described in Wechsler et al. (2002). 
In their paper they show how halo accretion histories are 
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Figure 8. Accretion histories for halo/3, where we show for each 
panel a different cosmological model. The fit to the accretion his¬ 
tories of each realization (solid yellow line) is obtained from Eq. 7. 


Figure 9. Cumulative number of subhalos with more than 400 
particles as function of their mass for haloa, halo^, haloy and 
halo(5 at 2 ; = 0 for each cosmology. 


Table 2. Values for the formation epochs ac for haloo:, halo/3, 
haloy in all four cosmologies. 


ACDM EXP003 EXP008e3 EXP006 
haloa 0.194 0.181 0.183 0.142 

halo/3 0.261 0.229 0.230 0.177 

haloy 0.225 0.210 0.215 0.173 


cosmologies the lower concentrations are not the result of 
formation histories but the effect of the modified dynamics. 


4.2 Subhalos 

In this Section we study the subhalo abundance, their 
radial distribution and circular velocities. 


well described by an exponential form that depends on one 
parameter, the formation epoch which is defined as the 
expansion factor at which the logarithmic derivative of the 
mass evolution falls below a critical value S. Specifically, the 
fitting form is given by 


M(a) 

M(l) 


= exp 




(7) 


with S' = 2. We used the fitting function of Eq. 7 to com¬ 
pute the formation epochs for each of the three halos for all 
cosmologies. As an example, we show in Fig. 8 the fitting 
result for halo/3 and we summarize the formation epochs for 
all halos in Table 2. 

As pointed out in Wechsler et al. (2002), Dutton & 
Maccio (2014) and Ludlow et al. (2013), in a ACDM cos¬ 
mology an early formation epoch leads to higher concen¬ 
trations. The same happens for dynamical dark energy cos¬ 
mologies, e.g. Klypin et al. (2003); Dolag et al. (2004). In¬ 
terestingly, in coupled dark energy cosmology this behavior 
is not preserved. Despite the fact that a stronger coupling 
brings to an earlier halo formation epoch, halo concentra¬ 
tions decrease when the coupling is increased, as clearly vis¬ 
ible from the EXP006 realizations. This phenomenon has 
also been showed in Baldi et al. (2010) and Baldi (2011), 
where an analysis of the modifications due to the coupling 
has been extensively carried out. The result is that the fric¬ 
tion term in Eq. 4 is responsible for making the halo ex¬ 
pand by altering its virial equilibrium through the injection 
of kinetic energy in the system, which in turns lowers the 
concentration (see Baldi 2011). This shows how in coupled 


4 . 2.1 Abundance 

The lower number of substructures present in EXP006 
halos compared to ACDM can be recognized in Fig. 3. In 
Fig. 9 we show the subhalo mass function, where only sub¬ 
halos that lie within R 200 and that have more than 400 parti¬ 
cles are considered. The total number of subhalos in EXP006 
realizations is always from 50 to 75% lower than in the re¬ 
spective ACDM cases, while differences between EXP003 
and EXP008e3 and ACDM are much less evident (~ 10%). 
Thus, the missing satellites problem (Klypin et al. 1999; 
Moore et al. 1999) can be progressively alleviated when in¬ 
creasing the coupling. Note that the differences in the sub¬ 
halos minimum mass among haloa, halo/3 and haloy, halo(5 
are due to the different resolutions used (see Section 3.2). 


4 . 2.2 Radial Distribution 

Fig. 10 shows the cumulative distribution of subhalos 
as a function of the distance from the main halo center nor¬ 
malized to the total number of subhalos within R 2 oo- All 
halos in all cosmologies show non-significant differences in 
the cumulative radial distribution. In order to better under¬ 
stand the distribution of subhalos, in Fig. 11 we show the 
differential distribution in a sphere of constant radius for all 
cosmologies, the radii are 350 kpc for haloa and halo/3, 250 
kpc for haloy and 200 kpc halo(5. The number of bins is kept 
the same for each halo in all cosmologies and the vertical 
lines show the virial radii. The distributions show a clear 
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Figure 10. Cumulative number of subhalos with more than 400 
particles as function of distance from the main halo center for 
haloa, halo/3, halo 7 and halo(5 at 2 ; = 0 for each cosmology. 




Figure 11. Differential number of subhalos in R 200 with more 
than 400 particles as function of distance from main halo center 
for haloo:, halo/3, haloy and halo(5 at 2 : = 0 for each cosmology. The 
vertical dashed lines mark R 200 for each halo in each cosmology. 
For a given halo, the binning is kept constant for all cosmologies. 


decrease of the number of subhalos in EXP006 halos com¬ 
pared to their respective ACDM cases, while for EXP003 
and EXP008e3 cosmologies differences are not so evident. 

As pointed out in Section 4.1.3, coupled cosmologies de¬ 
crease halo concentrations thanks to the presence of a fric¬ 
tion term in the equation for the evolution of density pertur¬ 
bations (Eq. 4), despite the earlier halo formation epochs. 
We claim that the same effect can be responsible for the 
lower number of subhalos compared to ACDM. Thanks to 
the extra friction term and to subhalo lower concentrations, 
subhalos that are falling into the main halo potential well are 
more heavily stripped, thus less halos with more than 400 
particles survive. If this claim is correct, we should be able to 
find a difference in the subhalos number distribution when 
we reach distances from the main halo center that are bigger 
than the radius from which the gravitational influence of the 


Figure 12. Number of subhalos in 3 R 200 with more than 400 
particles as function of distance from main halo center for haloa, 
halo^, haloy and halo(5 at 2 : = 0 for each cosmology. The ver¬ 
tical dashed lines mark R 200 and 2 R 200 for each halo in each 
cosmology. For a given halo, the binning is kept constant for all 
cosmologies. 


host halo is felt. In Fig. 12 we show the differential radial 
distribution of the number of subhalos out to about three 
times the virial radius of each halo. For the sake of clarity, 
we choose to show only the most extreme cases, EXP006, 
and ACDM for all four halos. The dotted lines represent one 
and two times R 200 for each halo in each cosmology. What 
we would like to stress, is that there seem to be a decrease 
in the number of subhalos in halos living in EXP006 cosmol¬ 
ogy compared to their ACDM realizations only within the 
gravitational influence of the main halo. Between 1.5 and 2 
R 200 this behavior inverts and halos living in the strongly 
coupled cosmology seem to have a larger or at least a compa¬ 
rable number of subhalos with respect to their ACDM cases. 
Thus, we ascribe the presence of a lower subhalos number 
to a massive stripping effect rather than EXP006 producing 
intrinsically a lower number of subhalos. We would like to 
stress on the fact that lowering the number of subhalos can 
also be achieved by warm dark matter cosmologies (e.g. An- 
derhalden et al. 2013), but the fundamental difference lies 
on the fact that those subhalos in warm dark matter cos¬ 
mologies were never formed, while in coupled dark energy 
cosmologies subhalos do form but they are heavily stripped. 


4-2.3 Circular Velocities 

Boylan-Kolchin et al. (2011) first showed that N-body 
simulations of a Milky-Way size halo predict a significant 
number of subhalos with circular velocities higher than the 
circular velocities that we measure for the brightest satellites 
of the Milky Way, which is surprising since these massive 
subhalos should not fail in producing stars. 

The discrepancy between ACDM prediction and obser¬ 
vations can be alleviated in multiple ways, starting from 
baryonic processes. Brooks & Zolotov (2014) suggest that 
baryonic feedback processes could be responsible for a dark 
matter redistribution, with the result of decreasing the cen¬ 
tral densities of the most massive subhalos. Rashkov et al. 
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(2012) point out that the possibility of star formation being 
stochastic below a certain mass would justify the Milky Way 
having massive dark satellites; furthermore, they highlight 
the fact that the tension between the Via Lactea II simula¬ 
tion and observations is only a factor of two in mass, which 
suggests that the uncertainty on the Milky Way virial mass 
could be a viable way out from the tension (Vera-Ciro et ah 
2013; Kennedy et ah 2014). Purcell V Zentner (2012) showed 
that there exists a significant variation in subhalo properties 
even when the host halos have the same virial mass. 

Last but not least, the discrepancy can be alleviated by 
appealing to non-ACDM cosmologies. The cases for warm, 
mixed (cold and warm) and self-interacting dark matter are 
considered respectively in Lovell et al. (2012), Anderhalden 
et al. (2012, 2013), Vogelsberger et al. (2012). In all cases 
they find that subhalos are less concentrated due to their 
late formation time, suggesting that alternative cosmologies 
can contribute to alleviate the tension between predictions 
and observations. 

In Fig. 13 we show the rotation curves for the twelve 
most massive subhalos at the moment of infall. We used 
the correlation between orbital energy and subhalo mass 
loss found in Anderhalden et al. (2013) to determine the 
subhalos ranking. Each row illustrates the twelve subhalos 
rotation curves for a given main halo in all considered cos¬ 
mologies, from the top down we show haloa, halo/3, haloy, 
halo(5. In yellow we show the observed values for i’circ(^i/ 2 ) 
for the brightest dwarf galaxies orbiting around the Milky 
Way, data are taken from Anderhalden et al. (2013) and ref¬ 
erences therein. Despite haloa and halo/3 having compara¬ 
ble masses, the tension between simulated curves and mea¬ 
sured points in the ACDM case is more evident in halo/3, 
supporting the fact that subhalo properties can vary even 
when host halos have the same virial mass (Purcell & Zent¬ 
ner 2012). The tension is alleviated in the case of haloy 
and even more halo(5, given their lower masses (Vera-Ciro 
et al. 2013). Overall, when looking at all halos in EXP003 
and EXP008e3 cosmologies these do not show significant 
improvement compared to their ACDM realizations in de¬ 
creasing the inner densities of subhalos. On the other hand, 
in the case of EXP006 cosmology, all four halos show such a 
dramatic decrease in subhalos rotational velocity peaks that 
rotation curves become incompatible with measured values. 
The dramatic decrease was to be expected given the choice 
of a large coupling parameter for EXP006 cosmology, but 
nonetheless it is useful to understand the effects of the cou¬ 
pling. 

4.3 Zooming-in on a dwarf halo 

To better explore the effects of the coupling at high res¬ 
olutions, we simulated a dwarf galaxy halo, haloe. We chose 
an isolated halo (no structures with comparable mass within 
four of its virial radii) and, given the results of Section 4.2, 
we only focused on the two most distant cosmological cases, 
ACDM and EXP006 cosmology. The virial masses are re¬ 
spectively 4xlO^M0 and dxlO^M©, with a mass resolution 
of 1.3 xlO^M©. We show in Eig. 14 the density maps for 
haloe in both cosmologies, ACDM on the upper panel. It is 
visible how the number of substructures decreases in the case 
with coupling. In Eig. 15 we show density profiles and rota¬ 
tion curves for haloe in both cosmologies. The effect of the 
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Figure 14. Density maps for haloe in ACDM (upper panel) and 
in EXP006 cosmology (lower panel). The side of each projection 
is 2 xR 2 oo- 




Figure 15. Density profiles and rotation curves for haloe in 
ACDM (black) and in EXP006 (red) cosmology. The inner min¬ 
imum radius is three times the softening length, while the most 
outer radius is four times R 200 respectively. The vertical dashed 
lines represent R 200 for each cosmology. 
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Figure 13. Rotation curves of the most massive subhalos at the moment of infall for each halo in each cosmology. From the top row down 
we show haloa, halo/3, haloy, halo(5, from left to right we show ACDM (black), EXP003 (cyan), EXP008e3 (blues), EXP006 (red). We 
estimate the subhalo mass ranking at the moment of infall using the correlation between orbital energy and subhalo mass loss found in 
Anderhalden et al. (2013). The yellow points are the observed values for r’circ(^i/2) for the brightest dwarf galaxies orbiting around the 
Milky Way. Data are taken from Anderhalden et al. (2013) and references therein. The masses of each main halo realizations is written on 
each panel. 


coupling is very evident in lowering the concentration and 
flattening the rotation curve. The values for the halo con¬ 
centrations are c = 15.2 and c = 6.5 for ACDM and EXP006 
cosmology respectively. Although the density profile in the 
coupled dark energy case is less concentrated, it is still cuspy, 
showing that in coupled cosmologies, as in ACDM, we are 
not able to produce a dark matter only cored density pro¬ 
file. The inconsistency with observation thus still persists, 
given the observational evidence that supports cored den¬ 
sity profiles for the satellites of the Milky Way (Walker & 
Peharrubia 2011; Amorisco & Evans 2012; Amorisco et al. 
2013). Interesting to note, by constructing a model in which 
both warm and cold dark matter are present and only the 
cold component is coupled to dark energy, a very high value 
(/3c ~ 10) for the coupling constant is favored (Bonometto 


et al. 2015) and simulated dark matter only dwarf halos show 
a cored density profile (Maccio et al. 2015). Fig. 16 shows 
accretion histories, the ratio of M 200 to its value today is 
plotted as function of scale factor. As in Section 4.1.3, we 
calculated the formation epochs as in Wechsler et al. (2002) 
and obtained Uc = 0.191 (ACDM), Uc = 0.146 (EXP006), 
confirming the finding that coupled dark energy models have 
earlier formation times, also for less massive halos. 


5 CONCLUSIONS 

We have performed the first study in coupled dark en¬ 
ergy cosmologies on high resolution simulations on galac¬ 
tic scales, with the aim to study the effects of the coupling 
between dark energy and dark matter on these scales, so 
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a = l/(l +z) 

Figure 16. Accretion histories for haloe in ACDM (black) and in 
EXP006 (red) cosmology. We show the ratio of M 200 at a given 
time and M 200 today as a function of the scale factor. 

far neglected in previous studies. We chose to investigate 
two viable models, one with constant coupling and one with 
varying coupling with redshift; we also chose a third case 
where the constant coupling value has been pushed beyond 
observational constraints to better investigate its effects. We 
then selected three Milky-Way size halos, a GxlO^^M© halo 
and a dwarf halo bxlO^M©, and studied their properties in 
a ACDM reference model and in the coupled cosmologies, 
resolving each halo with ~ 10® particles. 

We computed concentrations and formation epochs for 
all halos and we find that, despite the earlier formation 
epochs of the coupled cosmologies halos, these have lower 
concentrations. In a ACDM or a dynamical dark energy sce¬ 
nario, earlier formation epochs would imply higher concen¬ 
trations, but in the coupled dark energy case the reason 
for lower concentrations is not related to formation histo¬ 
ries, but rather to the modified dynamics. In fact, the equa¬ 
tion for the linear evolution of density perturbations (Eq. 4) 
shows the friction term —an extra term compared to 
the ACDM case that redistributes the dark matter particles 
and lowers the central densities, despite the earlier forma¬ 
tion times (see Baldi 2011). We find that this behavior is 
reproduced for all mass scales that we have investigated. 

In particular, subhalos can also be significantly less con¬ 
centrated. When falling towards their host, they are more 
heavily stripped once they start feeling the gravitational 
influence of the host halo. This translates into decreasing 
the number of subhalos compared to the ACDM realiza¬ 
tion and, additionally, subhalos are themselves less massive 
and less concentrated. For these reasons, coupled cosmolo¬ 
gies can be helpful in alleviating satellite-scales inconsisten¬ 
cies of ACDM. On the other hand, we find that in order to 
try to solve these issues with the coupling alone, one needs 
to use an extreme value for the coupling constant that is 
ruled out by observational constraints. In fact, only in the 
case with the highest coupling value the number of subha¬ 
los is significantly reduced (up to 75% less subhalos) than 
in the respective ACDM cases, while for the viable coupling 
cosmologies the decrease is much less significant (10% less 
subhalos). Moreover, we find that the distribution of the 
subhalos inside the main halo virial radius does not vary 
significantly among cosmologies. Lastly, less concentrated 
coupled cosmologies subhalos can in principle be useful to 
reconcile the inconsistency between the observed properties 
of the Milky Way dwarf galaxies and ACDM simulations pre¬ 


dictions, but once more a high enough value for the coupling 
must be assumed. Interestingly, allowing the introduction of 
massive neutrinos does alleviate the constraints on the cou¬ 
pling (see e.g. La Vacca et al. 2009), leaving coupled dark 
energy models dynamics on sub-galactic scales an interesting 
option. 

Overall coupled dark energy models can be as effective 
as ACDM in reproducing observations on sub-galactic scales 
and, for specific choices of the coupling, they can improve 
the agreement between predicted and observed properties. 
Hence, coupled models would need to be further investi¬ 
gated, possibly taking into account the effects of baryons 
at sub-galactic scales, which, as already shown in dynami¬ 
cal dark energy models (Penzo et al. 2014), are expected to 
amplify differences observed in the dark matter only case. 
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